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Capillary rise is one of the most well-known and vivid illustrations of capillarity; however, there 
is no solution as yet except certain segmental solutions for asymptotic regimes. This paper used the 
singularity-free Bush equation and successfully obtained its Taylor’s series solution. The solution 
revealed that capillary rise dynamics is mainly controlled by the Bond number and the Galileo 
number, while the Bond number is a key parameter within the solution. Due to the poor rate of 
convergence of the series solution, an approximate analytic capillary rise h(t) was proposed, which 


has been verified numerically. 
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INTRODUCTION 


Capillary rise is one of the most well-known and vivid 
illustrations of capillarity. Knowledge of capillarity laws 
is important in oil recovery, civil engineering, dyeing of 
textile fabrics, ink printing, and a variety of other field- 
s. It is the capillarity that brings water to the upper 
layer of soils, drives sap in plants, or lays the basis for 
the operation of pens [1-25, 27-33, 35-38]. The capillary 
action has been studied by a number of masters in sci- 
ence. Leonardo da Vinci (1452-1519) recorded the effect 
in his notes and proposed that mountain streams may re- 
sult from capillary rise through a fine network of cracks. 
Jacques Rohault (1620-1675) erroneously suggested that 
capillary rise occurs owing to suppression of air circula- 
tion in the narrow tube, which creates a vacuum. Geo- 
vanni Borelli (1608-1675) demonstrated experimentally 
that h ~ 1/a. Geminiano Montanari (1633-87) attribut- 
ed circulation in plants to capillary rise. Francis Hauks- 
bee (1700s) conducted an extensive series of capillary rise 
experiments, as reported by Newton in his Opticks, but 
was left unattributed. James Jurin (1684-1750) indepen- 
dently confirmed h ~ 1/a; hence Jurin’s Law. Albert 
Einstein’s first paper that was submitted to Annalen der 
Physik in 1900, was on capillarity [1]. 

Washburn [3] developed the equation to describe the 
rate of the penetration of liquids into small cylindri- 
cal capillaries, based on the laws of hydraulics. He 
assumed that the travelling distance in the initial tur- 
bulence period was negligible, while the Poiseuille re- 
gion covered practically the entire flow. According to 
Poiseuille’s law, when neglecting the air resistance, the 
rate of flow through a cylindrical tube assumed the fol- 
lowing form: Q = a (at + Lsa?), where a is the radius 
of the tube, h is the capillary penetration height, and 
Ls, is the slip length. Q is the rate of flow and for a 
cylindrical tube it is given by Q = ma?44. The total 
effective pressure AP consists of the external pressure 
P, the hydrostatic pressure P, and the capillary pres- 
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sure, given by the Young-Laplace equation AP = 450955, 


Combining the above relations and neglecting the s- 
lip length gives the so-called Lucas-Washburn equation: 
ap. = (Pa + Py + Bacos), If there is no external pres- 
sure, a commonly used form of Lucas-Washburn equation 
is Suh st + pgha? siny = 2aocosh, where w is the incline 
angle of the tube with respect to the horizontal surface. 
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FIG. 1: The incline capillary rise 


When the gravitational force is negligible, a well- 
known Washburn’s law can be obtained with the initial 
condition h(0) =0: h= rae which predicts burst- 
like behavior with the velocity being infinite at zero time. 
This singularity of the solution highlights a deep inconsis- 
tency of the above equation [20]. Rideal [4] also obtained 
this equation through a rather simple derivation. 

Taking into account the momentum of the liquid in the 
tube and end-effect drag on the fluid entering the tube, 
Brittin [6] derived a more rigorous formulations of the 
Lucas-Washburn equation, as follows: hae + $(@)? + 
seo uh +gh= Zo cos, where p is density, u is the 
viscosity, ø is the liquid-air surface tension o = yzv, 0 is 
the wetting angle of the liquid, h is the height of capillary 
rise, a is the capillary radius, and g is the acceleration of 
gravity. This equation assumes the Poiseuille flow profile 
throughout the capillary. 

The most popular equation of capillary rise dynamics 
is noe + ($) + Eruh + gh = Zo cos0. Although 


dt pa2 
this equation has been widely used in the literature, it 
has a singularity of t = 0, namely h = dh — œ as 


t — 0, which would lead to an ill-posed problem. This 
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equation cannot even deal with the natural initial con- 
ditions h(0) = h(0) = 0. A formal remedy is to take 
h(0) = \/20 cos 6/(pa) [5], neglecting such a logical draw- 
back as the acceleration of the liquid front at zero time is 
infinite, hence hence more criticism of this equation can 
be found in [8, 20, 31]. 


This singularity problem was firstly pointed out by 
Szekely et al.[8], who successfully removed the singulari- 
ty problem by composing the correct energy balance for 
the entry flow: (h+ Za) $} + 1.225(4)? + Sh + gh = 
me cos 6. Xiao et al.[28] considered the entrance pressure 
loss effects by using Dreyer’s model [13] and the dynam- 
ic contact angle effect by Newman’s model [7] to further 
modify Lucas-Washburn’s equation: (h + 1.028a) #2 + 
0.958( 42)? + Ži (h +0.25a) $ + gh = 22 cos. Bush [36] 
a a similar equation with a different coefficient of 
(a) 


Despite researching capillary rise dynamics for cen- 
turies, no solution has yet been found. One should fis- 
r obtain information for asymptotic regimes. Matching 
different solutions between the regimes has never been 
solved and, therefore remains an unsolved problem. It 
would be natural to think that one cannot satisfy those 
asymptotic solutions, and should hence seek a complete 
solution for the time domain. The following section fo- 
cuses on Bush’s model, including its series solution and 
approximate analytic solution. 


To make the study self-contained, it is organised as 
follows: following the introduction, the next section 
presents Bush’s formulated equation. Its series solution 
will then be obtained, and an approximate analytic solu- 
tion will be proposed by using the initial conditions and 
the series solution. This will be followed by the numerical 
validation and as well as discussion thereof, and finally 
the study will be concluded. 


BUSH’S FORMULATIONS OF CAPILLARY RISE 
DYNAMICS 


In physics the capillary dynamics process is a struggle 
between surface tension with the combining effect of wall 
resistance and gravity. In the initial stage the surface 
tension plays a leading role, and later the combined effect 
of wall resistance and gravity takes over. 


Capillary rise dynamics is governed by a combination 
of gravity, viscosity, fluid inertia and dynamics pressure. 
In respect of continuum mechanics, the dynamics pro- 
cess must satisfy the conservation law of momentum: 
4 m(t) 2] = F + [,pv-vu-ndS, where the second 
term on the right-hand side is the total momentum flux 


na2ph? = rh, hence the force balance on the column 


capillary may be expressed as 


(pra?h+ ma )h= 2naocos 
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inertia added mass capillary force 
1. 
— pra?hg — ra? zph? — 2Qnah-t . (1) 
Bme —— 
weight =y viscous force 


dynamical pressure 


We can estimate the added mass Mma from the change in 
kinetic energy as column rise from h to h + Ah, namely 
AE; = 4A(mU?), where m = me + Mmo + Meo (mass 
in the column, in the spherical cap, and all other mass, 
respectively). In the column, me = ra?hp, u = U; in the 
spherical cap mo = 27 a3 p; and in the outer region radical 
flow extends to oo, nut u(r) decays. Volume preservation 
requires: ta?U = 27a?u,(a) leads to u,(a) = U/2. Con- 
tinuity gives 27a7u,(a) = 27r?u,(r) and leads to u,(r) = 
a ur(a) = 4U. Thus, the kinetic energy in the far field 
is moU? = $f u(r}dm = $f u,(r)*p2ar2dr, 
which leads to mao = $ pra. Hence the change in ki- 
netic energy will be 


1 1 
OER = zA Mme + mo + Meo )U? + gm2U Au 
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The added mass is estimated as Ma = Mo +t Moo = Tra? p. 


For a typical capillary rise, a = 1074 ~ 5 x 1074m, u = 


ba 
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FIG. 2: The dynamics of capillary rise 


1.5 x 10°88, p = 750kgm~* and top velocity can be 
at V = 0.04 ~ 0.778m/s, hence the typical capillary 
rise Reynolds number Re = e = 10 ~ 195. At such 
low Reynolds number, if the tube is smooth one the flow 
in the tube is fully developed Poiseuille (laminar) flow, 
hence the viscous force is T, = -Hh 

Considering the above relations, the Bush equation of 
capillary rise dynamics was realised [36]: 


7 @h 1 dh 8u , dh 
h4 2 
a ae gage) 7 ona as 
with initial height and velocity boundary conditions of 
h(0) = 0, h(0) = 0. 


2 
+gh= £T cos 0. (3) 
pa 


Although no solution of Eq.(3) was obtained, some seg- 
mental solutions (as shown in Figure 3) for asymptotic 
regimes were proposed by a reduced form of Eq.(3). 
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FIG. 3: Asymptotic regimes of capillary rise dynamics [36] 


(1) Initial regime: h ~ 0, h ~ 0, infers h(t) = 22295, 
hence h(t) = ee (2) The quasi-steady stage 
(Lucas-Washburn equation), where the capillary force 
was compensated by gravity and the viscous drag, the 
asymptotic solution corresponding to the short-time lim- 
it (t > 0) is given by h(t) = 
limit (Quéré equation), the capillary rise was given by 
h= 20 cos 0 . 

pa 
ary level, ans was established by the balance of gravity 
and capillarity, hence H = are 

However, those segmental solutions were not matched 
into a single solution that was valid for the whole time 
domain. 


ag cos Ot 


n (3) Low-viscosity 


(4) The capillary liquid rose to a station- 


SOLUTION OF BUSH’S EQUATION OF 
CAPILLARY RISE DYNAMICS 


Eq. (3) is a nonlinear ordinary differential equation, and 
until now no closed form solution has been found. Fortu- 
nately, it has Taylor’s series solution h(t) = Xp o ent”, 


which is as follows 


_ 6gcosð > 30 cos0(Tgpa? + 240 cos 8) 


h(t) = t 
w Tpa? 343p7a® 
864107 cos? 65 é 
— ——.—1t O(t 4 
1715p3a™ PI (4) 


The correctness of this solution was verified by using the 
symbolic mathematics software, Maple. To understand 
the effect of surface tension, gravity and viscosity, Bush’s 
equation (3 was transformed into a dimensionless form, 
as follows: 


dz 7 Bo ďz ead all cos 6 dz, i 
z z =1, 
dr? 12cosdr? 2`dr Bo-Ga J dr 

(5) 


and the initial height and velocity boundary conditions 
were changed to z(0) = 0, 2(0) = 0. In which, z = 


3 


h ee 20 cos 0 = H 
mT = p, H = ga T= Vier the Bond number 


2 
o = ° represents the ratio of gravity and surface 


tension, and the Galileo number Ga = ga represents 
the ratio of gravity and viscosity. 

To express a fair quantitative understanding, we can 
rewrite the series solution in terms of Bo and Ga as fol- 
lows 


Qualitatively speaking, Eq.(6) reveals that the Bond 
number, Bo, plays a vital role in capillary rise dynam- 
ics, because it stands for the ratio of gravity and surface 
tension. In the beginning (within the diffusion time ay, 
the surface tension competes mainly with gravity, while 
after the diffusion time the viscous drag will come up a- 
gainst the surface tension. In other words, from the start 
of the capillary rise, the Bond number, Bo, enters the 
process, followed the Galileo number, Ga. 

The solution in Eq.(4) provides an important result, 
namely the non-zero initial acceleration h(0) = + coose 
which will later be used in the proposal of the approx- 
imate analytic solution. With the above solution in 
Eq.(4), there is no longer a need to divide the capillary 
dynamics process into several asymptotic regimes, since 
the solution in Eq.(4) is valid for the entire time domain. 
From Egqs.(6 and 4) we can see that the Bond number, 
Bo, is vital for the solution. Without the added mass ter- 
m Ta, the equations will have zero Bond number, which 
will lead to infinite height and velocity. This is the source 
of singularity. 

The academic value of the series solution enables one 
to realise the importance of the added mass, which can 
void the singularity. However, the rate of convergence of 
this solution is slow and useless for practical application. 
It would be natural to propose a simple solution for the 
capillary rise dynamics. 


APPROXIMATE ANALYTICAL SOLUTION 


Generally speaking, the entire dynamics process can 
be qualitatively described as follows: in an infinite reser- 
voir the capillary rise and the velocity are zero in the 
initial state. Due to the effect of the surface tension, 
the capillary liquid obtains initial acceleration (the ini- 
tial acceleration must never be zero), and begins to rise 
at a relatively uniform velocity, while the surface tension 
plays a dominant role in the ascending phase; however, 
as the capillary rises, wall frictions and gravity begin to 


work in an attempt to prevent the rise of the capillaries, 
and their joint action succeeds in decelerating the capil- 
laries at a point until the capillaries are finally stopped. 
Surface tension and wall resistance as well as gravity to 
achieve unity of opposites, and the capillary dynamics 
process are over, and attributed to calm. 

From a physical perspective of the capillary rising phe- 
nomena, the capillary liquid will rise to a stationary lev- 
el, established by the balance of gravity and capillarity, 
while the height and velocity will tends to h > H = 
20 cost and h — 0 at t + oo, respectively. To satisfy the 
conditions, we can propose that the capillary rise takes 
the following form 


h(t) = H[1 — fe), (7) 


where the exponent ( and function f(t) should be deter- 
mined. 

Substituting the Eq.(7) and h(t) = —Hf'(t)e~** + 
Hf (t)e~* into the initial conditions h(0) = 0, h(0) = 
0, while the initial condition for the function f(t) is given 
by 


f(0) = 1, f'(0) = £. (8) 


There are many possibilities for the construction of f(t), 
which can satisfy the condition in Eq.(8); for instance, 
f(t) = 1 + 6t+ att, in which the exponent d > 2 and a 
is a constant. Since the capillary rise dynamics process 
usually takes a millisecond, the lower order of time ¢ will 
dominate the process, and based on this understanding, 


we set the function f(t) = 1+ t+ at?. Hence, the 
capillary rise can be proposed as follows 

A(t) = H[1 — (1+ Bt + at? )e"**]. (9) 
Thus, we have capillary velocity h(t e H(aBt? — 3at + 


B2)te-Bt and arreleration h(t) = 
b6abt? + bat — B?)eP*. 

The capillary rise in Eq.(9) satisfies both initial con- 
ditions and infinity conditions. The unknown constants 
8 can be determined by the Taylor’s series solution in 


A(apt? + pt 


Eq.(4), which gives the acceleration h(0) = 2 coos? H, 

Hence, by equaling h(0 = H?, we obtain the exponent 
12 gcos@ _ 7 (10) 
7H pga? 


To find the parameter @ we can expand the capillary rise 
in Eq.(9) in the Taylor series, and obtain its coefficient 
of the fifth order team of t? as — 35H§?(6? + 15a), and 
let it equal the coefficient of t’ in Eq.(4), namely 


864 uo? cos? 0 
1715 


-ŻE + 15a) = , (11) 


paT 


which gives 


3 
Tuo cos 0 1 6g 
x , 12 

12p?a® 15 ( a 2) 


It is clear that a is crucial parameter, which links the 
capillary rise dynamics, to all constants. It is worth men- 


tioning an interesting aspect, namely that œ can be set 
22 2 
to zero by adjusting the radius as a = (EET, 


Finally, we obtain the capillary rise h(t) 


2 
E fı- ie (v2) ttal 


pga 
and the capillary velocity h(t) 


a one at? 89 _ 3at + 89 te (VT) 
dt pga 7a 7a 


dh(t) 
(14) 
as well as the capillary acceleration h(t). 

The above solutions in Eqs.(13) and (14) are valid for 
the entire time domain t € [0, o0), which to the author’s 
knowledge, has never been proposed in the lieterature. 

The solution in Eq.(13) reveals that the capillary rise 
h(t) is mainly controlled by the H = 725%? with the 
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decay rate of e TaJ where the radius a is the on- 
ly dominate parameter, and the smaller radius the fast 
decay. 


NUMERICAL VALIDATION 
Validation of the solution in Eq.(13) 


To validate the correctness of the solution in Eq.(13), 
we can compare the result obtained numerically by Zh- 
mud et al.[20]. For a comparative study, the data of the 
diethyl ether in glass is shown in the table below I. 


TABLE I: The diethyl ether in glass with arbitrary radius 
kg 


[22] o [53] 0 91%] p[2% 5] 
2.2. 10-7 1.67- 10-2 260 9.81 710 


In this case, for any radius a, we have the capillary 
rise H = 431 - 1076, hence the analytical capillary rise is 


given by 
4.31 - 107° 
h(t) = D 
2.9 3.821- 1071? 1.626, 3, — 20; 
[1+ ua z5 PETE Jte v=}. (15) 


If we set the radius a = 0.0005, we have the capillary rise 
as follows 


h(t) = 0.00862 [1 — (1 + 129.68¢ — 23112.54t?)e7 179-8] , 

(16) 
which is illustrated in the Figure 4. This figure shows 
that our solution agrees with [20]. 
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FIG. 4: Comparison study: a = 0.0005m and 0 = 26° 


The capillary rise velocity is given by 
h(t) = —25836¢(t + 0.06)(t — 0.08)e~129-®*, (17) 


which is illustrated in the Figure 5. The capillary rise 
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FIG. 5: Velocity profile for a = 0.0005m and 6 = 26° 


acceleration is given by 


h(t) = 3350496(t + 0.006) (t — 0.008) (t — 0.096)e~ 19-6°". 
(18) 


which is illustrated in the Figure 6. 


Influence of radius change 


By maintaining the wetting angle as 0 = 26° [20], and 
adjusting the radius of the glass in Eq.(15),reveals the 
graphs that are shown in Figure 7, which shows that the 
capillary rise is sensitive to the radius. The smaller the 
radius in which the liquid can travel, the further it goes. 
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FIG. 6: Acceleration profile for a = 0.0005m and 6 = 26° 


Therefore, the glass tube radius is a crucial parameter 
for capillary rise dynamics. 
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FIG. 7: The capillary rise of four different radius 


The radius also has a strong influence on the capillary 
velocity in Figure 8 and the acceleration, which can been 
seen in Figure 9. The smaller the radius, the faster the 
capillary moves. 


DISCUSSION AND CONCLUSION 


Piecing all the pictures together, we see a fair clear 
picture of the capillary rise dynamics as follows: the cap- 
illary rise onsets from an initial zero height and velocity, 
at t = 0, the surface tension, as the first driving force, 
supplies a kick-off acceleration of H82. This surface ten- 
sion steadily drives the capillary rising at an almost u- 
niform speed until it reaches its peak at a certain point; 
then the capillary speed is gradually reduced to zero after 
short oscillating owing to the combining effect of the wall 


capillary rise velocity 
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FIG. 9: The capillary rise acceleration for four different radius 


viscous friction and gravity. At the same time, the cap- 
illary acceleration decreases and reaches its lowest point 
before it increases to its final static state. 

Since the capillary rise generally has a low Reynolds 
number, the flow in the tube can be considered as a lam- 
inar flow. To avoid the singularity problem, the added 
mass must appear in the governing equation. From the 
dimensionless form of the Bush equation, we found that 
the capillary rise dynamics are mainly controlled by the 
Bond number and the Galileo number. We successfully 
obtained the Taylor series solution. Owing to the poor 
rate of convergence of this solution, we proposed a simple 
approximate analytic solution by using the Taylor series 
solution and initial conditions. We verified the proposed 
analytic solution by some numerical examples. 

Last, but not least, it may be worth pointing out 
that all previous solutions were from the non-oscillatory 
regime, however, when the liquid surface reaches h = H, 
the surface will oscillate before asymptotically ending. 


Within the oscillatory regime, the above solutions should 
be modified, which will be investigated by the Homotopy 
Analysis Method [26, 34] in the forthcoming paper. 
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